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ABSTRACT 

Aims. The EROS-2 project was designed to test the hypothesis that massive compact halo objects (the so-called "machos") could be 
a major component of the dark matter halo of the Milky Way galaxy. To this end, EROS-2 monitored over 6.7 years 33 x 10* stars in 
the Magellanic clouds for microlensing events caused by such objects. 

Methods. In this work, we use only a subsample of 7 x 10* bright stars spread over 84deg^ of the LMC and 9deg^ of the SMC. 
The strategy of using only bright stars helps to discriminate against background events due to variable stars and allows a simple 
determination of the effects of source confusion (blending). The use of a large solid angle makes the survey relatively insensitive to 
effects that could make the optical depth strongly direction dependent. 

Results. Using this sample of bright stars, only one candidate event was found, whereas ~ 39 events would have been expected if 
the Halo were entirely populated by objects of mass M ~ O.4M0. Combined with the results of EROS-1, this implies that the optical 
depth toward the Large Magellanic Cloud (LMC) due to such lenses is t < 0.36 x 10"' (95%CL), corresponding to a fraction of the 
halo mass of less than 8%. This optical depth is considerably less than that measured by the MACHO collaboration in the central 
region of the LMC. More generally, machos in the mass range 0.6 x IO^'Mq < M < 15Mq are ruled out as the primary occupants of 
the Milky Way Halo. 

Key words. Galaxy:halo - Cosmology: dark matter - Gravitational Lensing 
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1. Introduction 

Since the proposal (iPaczynskil 1 986t IPetroul 1981 *) that dark mat- 
ter in the form of faint compact objects (macho£]) could be 
found through gravitational microlensing, the EROS, MACHO, 
OGLE, MOA and SuperMACHO collaborations have monitored 
millions of stars in the Magellanic Clouds to search for mi- 
crolensing events. Such events would be due to a lensing object 
passing near the line of sight toward a background Magellanic 
star, causing a transient magnification of the star's primary im- 
age as well as creating a secondary image. Neither the image 
separation nor the image size are normally resolvable, so the 
only easily observable effect during an event is an apparent tran- 
sient amplification of the star's flux. The amplification is greater 
than a factor 1.34 if the line of sight to the star passes within the 
lens's Einstein ring of squared radius = 4GMDsx(l - x)/c^ 

' for "Massive Astrophysical Compact Halo Objects" ( lGriest|[r99 ll) , 
not to be confused with the "MACHO collaboration." 
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where is the distance to the star and xD^ is the distance to 
the lens of mass M. The optical depth for microlensing, i.e. the 
probability that at a given time a given star is amplified by more 
than a factor 1.34, is 
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where p is the mass density of lenses. For source stars in the 
Magellanic Clouds, the order of magnitude of t is fv^^Jc^ ~ fx 
10"^ where Viot ~ 220 km s"' is the rotation velocity of the Milky 
Way and / is the fraction of the halo mass that is comprised 
of lensing objects. The factor of proportionality between r and 
fv^^Jc^ depends on the structure of the Halo. The benchmark 
value is often taken to be that for a spherical isothermal halo of 
core radius 5 kpc, the so-called "S model" used by the MACHO 
collaboration jAlcock et al.ll2000bl: lGriestlll99lb . For the Lai-ge 
Magellanic Cloud (LMC) this gives 



: 4.7/ X 10" 
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For the Small Magellanic Cloud (SMC), the S model gives 
Tsmc ~ l-4Tiinc- For a flattened halo, one find s a smaller value, 
typically Tsmc ~ Time (ISackett and Gouldlll993h . 

Magellanic stars can also be lensed by non-halo stars, ei- 
ther in the Magellanic Clouds or in the Milky Way disk. Lensing 
b y disk stars i s expected to have an optical depth of order 10"** 
dAlcock et al.l i2000 b). The optical depth for lensing by "self- 
lensing", i.e. lensing by stars in the Clouds, is expected to range 
from ~ 5 X 10 in th e center of flie LMC bar to ~ 0.5 x 10 
at 3deg from the bar (iMancini et al.l 12004 ). For the SMC, the 
self-lensing optical depth is expected to be somewhat larger, 
~ 4 X 10~^ averaged over the central lOdeg^ dGrafi" & Gardiner! 
119991) . 

Microlensing events are characterized by a timescale fn giv- 
ing the time for the lens to travel a distance corresponding to its 
Einstein radius, - r^lv, where v, is the lens's transverse ve- 
locity relative to the line of sight. For high amplification events, 
2fE is the time over which the amplification is A > 1 .34. Since te 
is proportional to the square root of the lens mass M, the mean 
t-E will scale like M'^^. The S model has a 3-dimensional macho 
velocity dispersion of 270 km s ' and gives 
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Much excitement was generated by the MACHO collabora- 
tion's measurement of the LMC microlensing rate which sug- 
gested that a significant amount of the Milky Way's Halo is 
comprised of machos. Their latest analysis dAlcock et alJ2000bl) 
used 13/17 observed eventfl to measure an optical depth of 
Time/ 10"'^ = l-2!o3 istat) with an additional 20% to 30% sys- 
tematic error. This would correspond to a Halo fraction 0.08 < 
/ < 0.50 (95% CL). The mean of their events was 40 d corre- 
sponding to machos in the mass range 0.1 5Mo < M < 0.9Mo . 
On t he other hand the EROS collaboration (Lasserr e et al.ll2000t 
lAfonso et al. 200 3a) has placed only an upper limit on the halo 
fraction, / < 0.2 (95% CL) for objects in this mass range, rul- 
ing out a large part of the range of / favored by the MACHO 
collaboration. 

(Bennett (2005) ai-gued that the MACHO optical depth should 
be reduced to rime/ 10"^ = 1.0 + 0.3 in order to take into account 



^ 13 of the 17 events satisfied their so-called A criteria intended to 
identify high signal-to-noise events. The other 4 events, so-called B 
events, are selected by looser cuts. 



contamination by variable stars. This p aper made use of the ob- 
servation by the EROS collaboration dTisserandl l2004t) of fur- 
ther variability of one of the MACHO A candidates, indicating 
intrinsic stellar variability. The paper also noted that the spec- 
ti-um of the MACHO B candidate MACHO-LMC-22 indicated 
that the source is an active background galaxy, as reported in 
lAlcock et al.l d2001b ') where the event was eliminated from the 
sample for studying high-mass lenses. Using four MACHO A 
candidates whose microlensing nature was confirmed by preci- 
sion photometry and the one A candidate rejected as a variable 
star, Bennett (2005) performed a likelihood analysis to argue that 
1 1 ± 1 of the 13 A candidates are likely to be microlensing events, 
yielding the revised optical depth. 

Machos can also be searched for by monitoring M31 and 
looking for temporal variations of surface brightness consistent 
with a star in M31 being mic rolensed. Candidate events have 
been reported by the VATT dUglesich et al.1 12004|). W eCAPP 
(iR iffeser et all |2003|). POINT-AG APE dCalchi N ovati et aL 
20051), MEGA dde Jong et all l2006l) and Nainital (Joshi etal 
2005) collaborations. The POINT-AGAPE and MEGA col- 
laborations presented efficiency calculations allowing them 
to constrain the content of the M31 and Milky Way halos. 
The disagreement between these two collaborations parallels 
that between the MACHO and EROS collaborations with the 
AGAPE collaboration finding a halo fraction in the range 
0.2 </ < 0.9, while the MEGA collaboration finds a halo 
fraction / < 0.3. 

In this paper, we extend our previous analysis to find Time < 
0.36 X 10"^ (95% CL) for M ~ OAMq, corresponding to 
/ < 0.08. Unlike the previous EROS limit, this is significantly 
lower than the optical depth measured by the MACHO collabo- 
ration. Unlike all previous analyses, we use only a bright, well- 
measured subsample of the Magellanic stars, about 20% of the 
total. We believe that the use of this bright subsample gives more 
reliable limits on the optical depth than measurements using 
faint stars. There are two reasons for this. First, bright stars have 
well reconstructed light curves that permit discrimination of in- 
trinsically variable stars. Second, the use of bright stars makes it 
relatively simple to estimate so-called blending effects where re- 
constructed fluxes can receive contributions from more than one 
star, complicating the interpretation of events. 

EROS-2 is a second generation microlensing experiment. 
The first generation, EROS-1, consisted of two programs, both 
at the European Sou thern Observ atory (ESO) at La Silla, Chile. 
The first program (lAnsari et al.l il996a) used Schmidt photo- 
graphic plates to monitor a 27 deg^ region containing the LMC 
bar during the southern summer from October, 1990 through 
April, 1993. With a sampling frequency of up to one image per 
night, it was sensitive mostly to mac hos in the range 10 ^ ^Mp < 
M < IMq. The second program (iRenault et al.l ll998f) used 
a 0.4 deg^ CCD mosaic from December 1991 through March, 
1995 to monitor one field in the LMC bar and another in the 
SMC. With up to 40 images taken per night, this program was 
sensitive mostly to machos in the range 

The results of these two EROS-1 programs are summarized in 
IRenault et all dl997h . 

The second generation program described here, EROS-2, 
used the Marly 1 meter telescope at ESO, La Silla. The tele- 
scope was equipped with two 0.95 deg^ CCD mosaics to monitor 
93 deg^ in the Magellanic Clouds, 63 deg- in the Galactic Bulge, 
and 28 deg- in the spiral arms of the Milky Way. The observa- 
tions were performed between July 1996 and February 2003 (JD 
between 2,450,300 and 2,452,700). 
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Besides the Magellanic results presented here, EROS-2 
has also published me asurements o f the optical de pth to- 
ward the Galactic Bulge (lAfonso et al. 2003b; Hamadach e et alj 
I2OO6I) . The measured optic al depth is in agree ment with that 
me asured by MACH O tPopowski et alJ l2005l) . and OGLE- 



2 jSumi et alJ l2006l) and with the predictions of Galactic 



models ('Evans & Belokurov' '2002'; Bissa ntz & Gerhard I200I 
iHan & Go uld 2003; Wood and Mao 2005). 

EROS-1 and EROS-2 overlapped with the MACHO program 
that monitored -13.4 deg^ of the LMC from July, 1992 through 
January, 2000. Three other wide-angle microlensing searches are 
now in operation: MOA0 (since August, 1998), OGLE-^ (since 
June, 2001) and SuperMACHC0 (since October, 2001). 

In this article, we report on the analysis of the full EROS-2 
data set (July 1996 till February 2003) toward the Magellanic 
Clouds. Our previous analyses reporte d on 5 years of dat a 
and 5 million stars toward the SMC ( Afonso et all l2003al) . 
and 2 years of data and 17 million stars toward the LMC 
(iLasserre et al.ll2000l). Ari update of th e LMC analysis was re- 
ported in iMilsztain et alj dioOlh and iLasserr i mod), which 
dealt with 25 million stars and 3 years. The limits determined in 
the previous analyses are refined in the analysis prese nted here. 
More details about the present analysis can be found in lTisserandl 
(120041) . 

The plan of the article is as follows. In Section |2] we recall 
the basics of the EROS-2 setup, give the general characteristics 
of the data sample and describe the data reduction steps used 
to produce the light curves. Section [3] presents the Bright-Star 
Sample of stars to be used in the measurement of the optical 
depth. In Section |4] we describe the selection criteria used to 
choose the microlensing candidates in the Bright-Star Sample. 
(Selection of events in the full sample is described in Appendix 
IaI ) Section |5] presents the final sample of selected events from 
the Bright-Star Sample as well as events found by relaxing the 
selection criteria. In Section |6] we discuss the status of former 
EROS-1 and EROS-2 microlensing candidates as well as those 
of the MACHO collaboration. In Section]?] we describe the com- 
putation of the EROS-2 detection efficiency. Section [8] presents 
the limit on the optical depth and on the abundance of machos in 
the Galactic halo by combining all EROS-1 and EROS-2 data. 
We conclude in Section |9] with a discussion of th e significance 
of the h mit, and a comparison with the results of lAlcock et alJ 
(I2000bh . 



2. Experimental setup, observations and data 
reduction 

The EROS-2 Marly telescope, camera, telescope op e ratiori s 
and data reduction are describ ed in iBauer et a U (Il997h . 
iPalanque-Delabrouille et alj (Il998l) . and references therein. Here 
we give only general information, and details or modifications 
that are specific to the present analysis. 

2.1. The setup and the data 

The Marly telescope is a one meter diameter Ritchey-Chretien 
(f - 5.14m), equipped with two wide angle CCD cameras. 
Each camera is a mosaic of 8 CCDs, 2 along right ascension 
and 4 along declination. Each CCD has 2048 x 2048 pixels 
of 15 X 15 jum^ size, corresponding to 0.602 x 0.602 arcsec^. 



^ http://www.phys.canterbury.ac.nz/moa/ 
* http://bulge.astro.princeton.edu/~ogle/ 
^ http://www.ctio.noao.edu/~supermacho/ 
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Fig. 1. Map of the EROS-2 LMC and SMC fields in equatorial 
coordinates. A total of 88 LMC and 10 SMC fields were moni- 
tored. The first number in each field is the field number and the 
second is the number of bright stars (as defined in Section O in 
the field in units of 10'*. The two shaded regions (the larger one 
centered on the LMC bar) are the 13.4 deg^ used by the MACHO 
collaboration to measure the optical depth CAlcock et al..2000bl) . 



Images were taken simultaneously in two wide passbands, so- 
called 7?eros centered close to the Ic standard band, and Bgros in- 
termediate between the standard V and R bands. While no results 
presented here depend on the photometric calibration, almost all 
of our fields could be calibrated using stars fro m the catalogs 
of th e Magellanic Clouds Photometric Survey (IZaritskv et al.l 
l2004i) . For 4.5 deg^, the calibrat ion was checked with the OGLE- 
II catalog (lUdalski et al.l2000bl) . To a precision of ~ 0.1 mag, the 
EROS magnitudes satisfy 



= Vj - o.4(yj - /c) . 



(4) 



The observed fields correspond to 0.95 deg^ each and they 
are twice as large in declination than in right ascension (1.38 x 
0.69 deg^). A total of 88 fields have been monitored toward 
the LMC and 10 toward the SMC. The positions of the fields 
are shown in Figure [T] The exposure times ranged from 180 s 
to 900 s. The fields have been sampled differently, according to 
their stellar density and distance from the optical centers of the 
LMC and the SMC. Broadly, there are three famihes of LMC 
fields : in the first 3 years of operation, 22 outer fields were sel- 
dom imaged, 22 fields were imaged about 100 times and the 
remaining 44 inner fields were imaged over 200 times. Later, 
from July 1999 on, all 88 LMC fields were imaged with a simi- 
lar sampling. The number of photometric measurements, per star 
and per band, ranges from 300 to 600 in the LMC fields. 

The ten SMC fields were imaged with a similar sampling, 
except in 2001 and 2002 when the inner six fields were imaged 
twice as often as the four outer ones. The number of photometric 
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measurements, per star and per band, ranges from 600 to 900 in 
these SMC fields. 

After rejection of bad images (11%) due to bad seeing, high 
sky background or instrumental problems, the average numbers 
of measurements used in the analysis, per star and per band, are 
about 430 and 780 in the LMC and SMC fields respectively. 

2.2. Template Images 

The production of light curves proceeded in three steps : tem- 
plate image construction, star catalog production from the tem- 
plates, and photometry of individual images to obtain the light 
curves. 

The template images were obtained by co-adding 15 images 
of each field; they were resampled so that the templates have 
twice as many pixels as the original images. The 15 images were 
chosen among the best ones available, i.e. with low sky back- 
ground, good seeing, and with a large number of stars (as esti- 
mated from a quick first look algorithm). All images susceptible 
to enter the template construction were checked for the absence 
of long tracks, caused by satellites, planes or meteors. For tech- 
nical reasons linked to computing and to PSF variation within 
one CCD, the CCD images were divided in four quadrants, such 
that there were in total 6272 template images (98 fields, 8 CCDs, 
4 quadrants, 2 passbands). In order to ensure relatively uniform 
zeropoints of the 6272 templates, we required that the first im- 
age used in template construction (to which the other 14 images 
were photometrically aligned) be registered within a short time 
interval with good and uniform sky conditions (23rd to 26th, 
November 2000). For each field, the same epochs were used in 
the construction of the templates in the two passbands. 

The stars were identified on the template images using a 
pseudo-image that we call a correlation image. Each pixel of 
this image contains the correlation coefficient of neighboring 
pixels of the template itself with a two-dimensional Gaussian 
PSF. Each group of pixels satisfying some threshold value on the 
pseudo-image was retained as a star in the catalog. In previous 
EROS-2 analyses, the thresholds were identical for all templates; 
this had led to over 20% failures in this cataloging step. The 
present analysis has chosen to progressively relax the thresholds 
when such failures occur; in this way, the cataloging step fail- 
ures were drastically reduced. (The number of identified stars on 
fields using relaxed thresholds was lower on average.) The star 
catalogs obtained from the template images in the two passbands 
were then merged. A star was retained only if it was detected 
in both. The overall efficiency of the template plus star catalog 
construction was excellent; only 24 CCD quadrants could not be 
processed out of a total of 3136. 

Examples of color-magnitude diagrams can be found in 
Figures [T] to [10] and lA.ll to IA.5I They are all characterized by 
a prominent group of clump giants and a main sequence whose 
relative strength varies from field to field. There are also stars 
that are redder than Magellanic red giants. Most of them are 
likely to be foreground stars in the disk of the Milky Way and 
thei r number is consis tent with the predictions of Galactic mod- 
els dRobin et alj|2003h . 



2.3. Light Curves 

Photometry was then performed on each image of a given quad- 
rant in turn with software specifically designed for crowded 
fields, PEIDA ( Photometrie et Etud e d'lmages Destinees a 
r Astrophysique) (jAnsari et al.ll996bl) . First, the image was geo- 
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Fig. 2. The photometric precision as a function of /?ero.s for the 
dense field lm009 (top) and the sparse field lm048 (bottom). 
Each point represents the rm.s. dispersion of the measurements 
for a single star after elimination of outliers. The vertical line 
shows the position of the bright-star magnitude cut (|5]l for the 
field, «eros = 18.23 for lm009 CCDs 0-3 and i^ems = 19.7 for 
lm048 CCDs 0-3. 



metrically aligned with the template. Then, imposing the star po- 
sition determined from the template image, PSF-fitting photom- 
etry was performed for all stars by a linear least-squares method 
involving the star and all neighbors closer than 1 1 pixels, plus a 
sky background. An estimate of the error on this flux measure- 
ment was computed that depends on the photon statistics and on 
the overall quality of the photometered image. Typical photo- 
metric precisions are shown in Figure |2] 

Before the analysis was started, we removed from the light 
curves measurements taken under far from normal conditions. 
This happened not infrequently, as the data taking policy was 
to work whenever possible, leaving to the analysis the task of 
rejecting these abnormal measurements. These were identified 
by extreme values of the sky background, seeing or absorption. 
In addition, images where the photometry failed for over 40% 
of cataloged stars were eliminated, as well as images for which 
more than 12% of the stars showed an excursion from their aver- 
age flux larger than three standard deviations. Depending on the 
CCD, the number of rejected images varied between 7 and 18%, 
with an average of 1 1 %Q 



* The largest single cause for rejection was the malfunction that af- 
fected 5 CCDs of the R„o<i passband camera starting in January 2002. 
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To reduce systematic errors in the photometry, each light 
curve was searched for significant Unear correlations between 
the measured flux and three observational variables, the seeing, 
the hour angle and the airmass. This was done independently in 
the two passbands. The measured fluxes were corrected linearly 
by requiring a vanishing correlation coefficient between the cor- 
rected fluxes and the given variable. The largest correlation was 
found with the seeing, in both passbands. For the bright stars 
considered in this work, the correction has only a small influence 
on the point-to-point dispersion of the lightcurves, reducing it on 
average by 6%. 

A potential problem with this correction is that an artificial 
correlation can be induced if a real flux variation happens to oc- 
cur during a period of poor seeing (for example). In this case, ap- 
plication of the correction would reduce the amplitude of the real 
flux variation but increase the point-to-point dispersion of the 
curve by making an incorrect flux correction. To guard against 
this possibility, we did not apply the correction if it increases the 
point-to-point dispersion of the light curves. 

2.4. The Full Sample of Stars 

A total of 58.4 million objects were found on the template im- 
ages of both passbands - 51.8 million in the LMC and 6.6 mil- 
lion in the SMC. We chose to ignore those stars for which the 
association between the objects detected separately in the two 
passbands was doubtful or ambiguous (8.4 million) and objects 
that are dimmer than about twice the typical size of sky back- 
ground fluctuations (14.9 million). We rejected light curves for 
which more than half of the photometry points are absent (1.0 
million). Finally, we did not consider stars whose photometry is 
unstable due to its environment. This includes stars close to a 
very bright field star (Vj < 10.1, probably in the Galactic disk) 
and stars close to a visible diffraction feature in the PSF of bright 
stars; these two categories contain respectively 1.1% and 0.8% 
of the remaining stars. 

We removed from consideration stars in field lm003, 
CCDO which has numerous events caused by light echos from 
SN 1987a. The echos generate arc-like images that appear to 
move a few arcsec per yeafl causing false variations of a star's 
flux as the arc passes through the star's position. 

After these cuts, 33.4 million objects remained, 29.2 million 
in the LMC and 4.2 million in the SMC. This constitutes the 
Full Sample of EROS-2 stars. In the next section we describe the 
selection of the Bright-Star Sample used for the measurement of 
optical depths. 

3. The Bright-Star Sample 

We have chosen to restrict our analysis to the Bright Sample of 
stars defined below. Besides the obvious advantage of ensuring 
a good photometric resolution, we do this to simplify the eval- 
uation of the number of expected events predicted by a model. 
In crowded fields, this evaluation is complicated by "blending", 
i.e. the fact that photometry of a given object can be influenced 
by more than one star. We shall see that these effects are rather 
small and simple to evaluate for the Bright-Star Sample (which 
should, strictly speaking, be called the Bright-Object Sample). 



We have thus chosen to reject all measurements taken in 2002 and 2003 
for these 5 CCDs. 

' Using 1000 pictures of this region co-added by groups of 10, we 
have produced a film that shows the motion of the light echos; it can be 
found at the URL ht tp : //eros . in2p3 . fr/EchoesSH1987a/| 
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Limiting the number of studied objects limits the sensitivity 
of the experiment so the magnitude cut must be a compromise 
between quantity and quality of objects. The efficiency calcula- 
tions of Section |7] for unblended sources indicated that keeping 
only the ~ 20% brightest stars reduces by only ~ 50% the num- 
ber of simulated microlensing events that pass our selection cri- 
teria. We therefore initially decided to accept a nominal ~ 50% 
loss of events by requiring /?eios < ^median, where /^median is the 
median /?eios for simulated unblended microlensing events pass- 
ing our event selection criteria, ^median ranges from ~ 18.2 in the 
densest fields to ~ 20.5 in the sparsest. However, for sparse fields 
far from the LMC bar, in order to have a reasonable object re- 
construction efficiency, a stricter cut was found to be necessary, 

/?ei-os < 19.7. 

The final sample of bright stars is therefore defined by 

16.0 < T^ei-os < ^max ^max = ""'"(^median, 19.7) . (5) 

The minimum magnitude /?ei-os = 16 was chosen to avoid the 
large number of variable stars brighter than this. 

The position of the magnitudes cut are shown in the color- 
magnitude diagrams of Figuresl7lto[T0land lA.ll to lA.5l Generally 
speaking, the cut includes clump giants but not the numer- 
ous main sequence stars seen far below the clump. Other than 
the small number of bright main sequence stars, we there- 
fore employ stars of colors and magnitudes similar to those 
used in Galactic Bulge measurements that use clump giants 
(Hamadache et al. 2006). 

It turns out that the cut (|5]l gives a rather uniform photo- 
metric precision for the Bright-Star Sample. Figure |3] shows the 
distribution of precision with a mean in both bands of ~ 7%. 
The precision in dense fields hardly differs from that in sparse 
fields and the precision in the LMC hardly differs from that in 
the SMC. 

The number of EROS-2 objects in the Bright-Star Sample 
defined by (|5]l is 6.05 million in the LMC and 0.90 million in the 
SMC. We must subtract from these numbers the expected num- 
ber of foreground stars in the Milky Way d isk. The Besangon 
model of the Milky Way jRobin et al.ll2003[) predicts that 0.56 
million (LMC) and 0.04 million (SMC) are foreground Milky 
Way stars, consistent with the number of stars we observe that 
are redder than Magellanic red giants. The number of Bright- 
Sample stars to be used for optical depth measurements must 
therefore be reduced by 9% (4%) for the LMC (SMC), giving 

N^bjecs^ 5 49x 10'' (LMC); (6) 
N„bjec,s =0.86x10'' (SMC). (7) 

We now show how we evaluate the effect of blending on the 
Bright-Star Sample. For a given lens model (defined by the spa- 
tial, velocity and mass distribution of lenses), the expected num- 
ber of detected events is 

A^ex = -^ry<e,) (8) 

where Tobs is the observing period, (t^) is the mean Einstein ra- 
dius crossing time given by the model, e.g. (|3]l, and r is the opti- 
cal depth. The sum is over all observed stars and the "efliciency", 
(e,), is the expected ratio between the number of events passing 
the selection cuts of Section |4] and the number of events with 
M < 1 occurring during the observation period. The efficiency 
must be averaged over the distribution given by the model. 
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Fig. 3. The flux precision, A(/>/0, for the Bright-Star Sample (2% 
of the sample are in the histograms). Each entry is the mean 
flux uncertainty for a star on its light curve. The top (bottom) 
histograms are for the 7?eios (Beros) bands. 



EROS photometry is performed on a set of "objects" found 
on a reference image. In crowded fields, the correspondence be- 
tween objects and stars is not straightforward. Generally speak- 
ing, most objects are dominated by one star, i.e. by the brightest 
star in the object's seeing disk. However, if a fainter star in the 
seeing disk is microlensed, EROS photometry will assign some 
of the extra flux to the object with the remainder assigned to sky 
background. It is thus convenient to rewrite ([SJ as 



(9) 



where the second sum is over the stars / near the object j. The 
efficiency, (e,j), is again the expected ratio (averaged over fe) 
between the number of microlensings of star / that pass selec- 
tion cuts and the number of microlensings of star / with u < I 
occurring during the observation period. Since we do not know 
the stars in the seeing disk of each object, we must estimate the 
sum statistically using knowledge of the luminosity function. We 
will see that sum is dominated by the brightest star, / = 1, with 
contributions from the second brightest, / - 2, of order 10%. 



If there were no blending, the efficiency could be estimated 
by modifying a sample of light curves with measured fluxes 
Fo(t) by 

^-(0 = ^^0(0-^^^^, u(tf^ul + [(t-to)/tE]\ (10) 
u Vm^ + 4 

where u(t) is the distance of the lens to the star's line-of-sight 
in units of the Einstein radius, mq is this distance at the time, 
to, of maximum amplification, and Je is the Einstein ring radius 
crossing time. The modified lightcurves can then be subjected 
to the selection criteria and the efficiency deduced by averaging 
over and uq. 

For blended events only a fraction a of an object's ffux is 
amplified. If the fraction is independent of amplification, the re- 
sulting light curve is 



F(t) = Fo(t) 



(I - a) + a- 



U-+2 



(11) 



u Vm2 +4_ 

The efficiency to see such events now depends on a: e/j - e,)(a,). 
Stars with small or,- require a small impact parameter mq to give 
a sufficiently large reconstructed amplification. Since our selec- 
tion criteria require a reconstructed amplification greater than 
1.34 (uq = 1 for a = 1), the primary effect of blending is to 
reduce the efficiency by a factor Uiia,), the impact parameter 
needed for star / to produce a reconstructed amplification of 1 .34. 
The efficiency is reduced by this factor because of the ff at a priori 
u distribution. Additionally, the event appears shorter by a factor 
M, further modifying the efficiency. 

A simple example is two superimposed stars of the same 
magnitude and color. In this case a\ - ai - 0.5 and u\ - 
U2 ~ 0.7. Before efficiency corrections, the event rate would be 
0.7 + 0.7 = 1.4 times the rate calculated ignoring blending. If 
the two stars have colors differing by A(B - R)eros = 0.7, the 
requirement that A,„a_^ > 1 .34 in both colors yields a rate that is 
increased by a factor 1.3 over the rate for a = 1. 

To statistically evaluate the distribution of the c,- for Bright- 
Sample objects, a spectrum of 8000 artificial stars was placed 
at random positions on real EROS images (lm009, lm019, 
lm034, lm048) to give artificial images (xm009, xm019, xm034, 
xm048). These fields were chosen as representative of the 
crowding variations over the LMC fields. The stars were given 
ffuxes according to randomly chosen microlensing events (one 
per artiffcial star). Reference xm images were created using these 
images and photometry performed in the same manner as for the 
normal (Im) images. 

After photometry, the light curves of xm objects falling near 
artiffcial stars were studied to ffnd values of the a,. In practice, 
one median value of a was calculated for each lightcurve us- 
ing points for which (A - l),„p„, > 0.5, though no nonlinearity 
was observed that made this a crucial point. Of the 8000 artiff- 
cial stars that were added to each of the four ffelds, most are not 
usable for various (understandable) reasons: star in a masked re- 
gion of the CCD, event during a period with no observations, 
or the reconstructed magnitude not satisfying the Bright Sample 
cuts. For the ffeld xm009, a total of 11 23 lightcurves were usable 
for determining a. Of these lightcurves, 982 concerned objects 
for which the artiffcial star is the primary component, while in 
the remaining 141 objects the artiffcial star falls underneath a 
preexisting bright object on the original image. We use the ffrst 
type of object to determine the distribution of ai, the value of a 
for the primary star associated with each object. 

Figure|4](top) shows the distribution of a[ for the dense ffeld 
lm009 and the sparse ffeld lm048. The other two ffelds give sim- 
ilar distributions. The distributions of ai are characterized by 
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Fig. 4. The effects of blending on artificial stars added to the dense field lm009 (left) and the sparse field lm048 (right). The top 
two panels show the distribution of ai using artificial stars that do not fall under a brighter pre-existing object. Stars further than 2 
arcsec from the nearest pre-existing object are concentrated in the peak at o-i ~ 1 while stars nearer to pre-existing objects form the 
tail at a I < 1. The middle two panels show, for objects within 2 arcsec of pre-existing objects, the correlation between ai and the 
difference in magnitude between the artificial star and the pre-existing object. Artificial objects with ai significantly below unity 
are associated with pre-existing objects of similar magnitude. The line shows the expected relation for two superimposed objects: 
ai - l./(l. + lO^"*^). The bottom two panels show the distribution of (01,02) with 02 calculated from (fT2b and the magnitude of 
the pre-existing object. The hne shows the expected relation for two superimposed objects: ai + 02 = 1. 



peaks at ffi ~ 1 due to artificial stars falling more than ~ 2 arcsec 
from any pre-existing Im object. This happens for about half the 
artificial stars in the densest field (lm009) so blending has little 
effect on about half the bright stars in dense fields. Stars falling 
on pre-existing stars yield the tail at ffi < 1. This can be seen in 
the middle plots that show, for artificial stars within 2 arcsec of 



pre-existing stars, the correlation between ai and the magnitude 
difference between artificial and pre-existing stars. The distribu- 
tion follows closely the expected form = 1/(1 + ./2//1) where 
/i is the flux of the artificial star and /2 is the flux of the preexist- 
ing star. For example, when ai ~ 0.5 the cloud of points passes 
near AR ~ 0. 
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The mean value of ai is ~ 0.86 equivalent to a mean value 
of Ml of 0.92. This means that the (efficiency uncorrected) rate 
of events with amplifications of the primary component greater 
than 1.34 is reduced by a factor 0.92 from the rate calculated 
assuming no blending. 

The loss of 8% of the events is compensated by extra events 
due to stars under the primary component. Since ~ 1/(1 + 
fi/fi), we can expect that 02 ~ 1/(1 + f\l fi) for superimposed 
objects. However, the superposition is not perfect so 02 must be 
a decreasing function of the separation between the object and 
star 2. The artificial events due to artificial stars falling under 
brighter pre-existing objects are well described by 



lmd48 



a2{d) 



fi 



A +/2 



exp[-(d/do) ] 



da = 1.9 arcsec (12) 



where fi is the flux of the major component, /2 is the flux of the 
minor (amplified) component and d is the separation between the 
xm object and the minor component. This formula can then be 
applied to the objects where the artificial star is the major com- 
ponent since the flux and position of the minor component (on 
the original image) is known. In the bottom panels of figure |4] 
we show ai vs. 02 for the fields lm009 and lm048. The anticorel- 
lation between a; and 02 means that the loss of events due to 
ai < 1 is partially compensated by events from the second star 
on an object by object basis. 

The compensation is seen in figure |5] histogramming ui and 
Ml + U2- The sum mi -n M2 gives the rate for event with amplifi- 
cations > 1.34 compared to that calculated assuming a = 1. In 
lm048, events due to the second star compensate for the loss of 
events on the first star while in the dense field lm009 there is 
a 10% overcompensation. As we will see in Section [T] the ef- 
ficiency to observe the blended events is degraded compared to 
unblended events, so the overall number of expected events is 
about 10% less than what one would calculate ignoring blend- 
ing. 

We have checked our calculations by studying public HST 
images of a small part of the densest EROS region, that of lm009 
ccd 3. A catalog of stars for this region was produced and images 
with EROS seeing were fabricated. As in the previous analysis, 
images with some amplifications were compared with reference 
images. Figure |5] shows the distribution of mi for bright objects 
on the convoluted images. Because of the very high density of 
stars. Ml is somewhat smaller than that calculated with the artifi- 
cial images: (mi) - 0.85 as opposed to (mi) = 0.92. This lower 
value of Ml is, as expected, compensated by higher values of U2- 
Figure|5]shows the histogram of mi + M2 indicating that the event 
rate is raised by 15% over the unblended rate. Including the loss 
of efficiency modifies this so that the event rate is only 5% higher 
than the unblended rate. 

4. Candidate selection from the Bright-Star Sample 

The present analysis aims at detecting luminosity excursions, 
due to microlensing, on otherwise constant light curves. The 
smallest reachable microlensing timescale, fe, is determined by 
the sampling of the fields; in practice, except for a few fields, 
the detection efficiency is very low below 2-3 days. The largest 
detectable timescale, about 800 days, is determined by the re- 
duction in the detection efficiency when the timescale becomes 
so large that the baseline flux is not seen during the observing 
period. 

The analysis is guided using the simulation of microlensing 
described in Section |7] which also serves to determine the effi- 
ciency of the selection procedure. 
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entries 982 
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entries 112 
<ul>= 0.844 
<ul+u2>= 1.175 
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Fig. 5. The distribution of mi (solid Une) and (mi + M2) dashed 
line for the sparse EROS field lm048 (top panel) the dense EROS 
field lm009 (middle panel) and the very dense HST field in 
lm009, ccd3 (bottom panel) as described in the text. The dis- 
tributions of Ml are all characterized by a peak near mi ~ 1 and a 
tail at M 1 < 1 . With increasing star density, the mean of m 1 -H M2 in- 
creases because of the increasing importance of secondary stars. 



In the present analysis, luminosity excursions are defined 
with respect to the baseline stellar flux; the first task is thus 
to determine this baseline. To that end, we order the mea- 
sured fluxes O; of the light curve by increasing values of O. For 
each of the N - I middle values of the flux intervals <]>,;„„d = 
0.5(<I), + <I>,+i), we count the number A^™./ of runs, i.e. groups of 
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cut 


T \/fr^ 




Bright-Star Sample 






Filter (C^i, L-Z, 


J Jo /O 


J lol 




6778 


847 


good (d^ 


2636 


266 


iQ \^\J} 






peak ^ ^ ' ) 


1624 


149 


significant A^-^ (C8) 


350 


32 


good4,„, (C9) 


152 


17 


BlueBumper(ClO) 


10 


4 


SN(5 <0.3) (CU) 


7 


4 


Mo < 1 (C12) 


1(SN) 


1 



Table 1. The number of light curves surviving the cuts CI to 
C12 as described in the text. 



100 




Fig. 6. The uq (mean of those for R^ms and Beros) distribution 
for LMC and SMC Hght curves surviving cuts C1-C7, C1-C8, 
and C1-C9. The positions of the 14 light curves surviving cuts 
CI -CIO are shown with the 4 supernovae as dashed lines. 



consecutive points on the light curve that are on the same side 
of <l>,,mid0 The baseline flux is defined as that value of <l>,,,„,y 
which maximizes N „,„,!■ This way of defining the baseline has 
a precision similar to that of a simple average of the fluxes, but 
it proves much more robust to aberrant measurements, and less 
biased for most variable stars. We have checked on the simu- 
lated microlensing light curves that there is no visible bias for 
timescales shorter than ~ 200 d. The bias is non-negligible 
above 600 d, but this is taken into account in the efficiency cal- 
culation. 

Next, we determine the point-to-point dispersion in the light 
curve cTptp, from the comparison of each measured flux with the 
linear interpolation of its two neighboring (in time) fluxes. This 
is done separately for the 7 seasons of data taking, as we have 
observed a progressive degradation of the photometric scatter in 
the last 3 seasons. 

Using the baseline flux and the photometric scatter, we then 
search for luminosity excursions defined as a group of consec- 
utive points with fluxes sufficiently far from the baseline flux. 
These should have at least 5 consecutive measurements more 
than 1 -ScTptp from the baseline. In practice, to increase the de- 



tection efficiency for short duration phenomena, we allow small 
"holes" within the excursion (series of points less than 1.5cr^,/p 
from the baseline) provided that each hole contains exactly one 
measurement. We call positive (negative) excursions those com- 
prised of fluxes higher (lower) than the baseline. 

The relative significance of each excursion is then estimated 
using the probability that it corresponds to a statistical fluctua- 
tion of a normal law. We use the variable LP^, the co-logarithm 
of this probability 

..„ = -|,„.(iErfc(^)). 

where Xj is the deviation from the baseline of the point taken at 
time f,, in units of its error cr,, and is the number of points 
within the excursion. The excursions are ranked by decreasing 
values of LP^ (decreasing significance). 

We now describe the 12 selection criteria, CI -CI 2, used to 
select microlensing candidates. CI requires at least one signifi- 
cant positive excursion: 



CI 



LPn.pos,i > 20 . 



The largest negative excursion should be much smaller than the 
largest positive one: 

C2 : LPN,pos.l/LPN,nes,l > 10 . 

There should be less than 10 excursions in total: 



C3 : 



excursions 



< 10 . 



Cuts CI, C2 and C3 are applied independently in the two 
passbands and a star is retained for further analysis only if it 
passes the three cuts in both bands. There are 33876 LMC light 
curves and 5787 SMC light curves that survive. The progression 
of the number of surviving light curves with the subsequent cuts 
is shown in Table [T] 

All light curves passing CI - 3 are fitted with a simple mi- 
crolensing curve ([Tol l, independently in the two passbands. The 
results of the fits, i.e. the fitted parameters fo, mq and and the 
values ofx^, are used in cuts C5-C10. 

Cuts C4 and C5 eliminate light curves with significant vari- 
ability outside the main positive excursion. We first require that 
the second positive excursion be much less significant than the 
first one: 



C4 : 



LP, 



N,pos, 1 



ILP. 



N,pos,2 



> 10 



Or, equivalently, Nrunj is the number of times the time-ordered light 
curve crosses the constant flux line at O,.,,,,^, plus one. 



To discriminate against low amplitude variations (either real or 
due to photometric problems) we require that the normalized 
xlase ^ independent light curve outside the main excur- 
sion (points with m > 2) be sufficiently good: 

C5 : < 20 . 

base 

This cut may seem loose as it would correspond to cuts at 20 
standard deviations, for a normal law. However, the errors we 
use to compute are only estimates of the actual eiTors. For 
isolated stars, the errors are overestimated. The values of the cuts 
have thus been tuned using the simulated light curves, that have 
the same uncertainties as the original data; the two cuts reject 
about 10% of the simulated microlensing light curves. 

Cuts C6 and C7 require that the main fluctuation occur at a 
time when the light curve is sufficiently well sampled. In order 
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to discriminate against stellar variations with very asymmetric 
peaks we thus require that the time fo be within the observing 
interval: 



C6 



Then + < to < T, 



end 



where Tt,eg and T^nd are the first and last date of EROS-2 obser- 
vations. We also require that there be more than 4 flux measure- 
ments within the main fluctuation (m < 2) 



C7 



N, 



peak 



>4 



There are 1624 LMC light curves and 149 SMC light curves 
that survive cuts C1-C7. Their distribution of mq and those of 
light curves satisfying subsequent cuts are shown in Figure|6] It 
is seen that most light curves show only low amplitude variations 
with Mo > 0.6. 

We next require that the microlensing fit to the light curves 
in both passbands be significantly better than the fit of a constant 
flux: 



C8 



Xct ' 



' ^ml 



1 



Xmll^dof yl2Ndof,peak 



>40 . 



where xlt Xmi the chi-squared values for the constant fit 
and the simple microlensing fits, and Ndof and Ndof,peak are the 
number of measured points in the full light curve and in the peak. 
The event is required to pass this cut in both passbands. This cut 
ehminates light curves that have low amplifications, poor sam- 
pling, and/or poor photometric errors. 

We then require a reasonably good fit to the microlensing 
curve within the peak: 



C9 : 



X\eak ^dof,peak 



^J2Ndof,peak 



< 10. 



There are 152 LMC and 17 SMC light curves that pass cuts 
C1-C9. As can be seen in Figure |6] almost all of these light 
curves have mq > 0.8. They are mostly bright main-sequence 
stars that display low-amplitude variations that are diflicult to 
distinguish from microlensing, though they often show ~ 20% 
more variation in the red than in the blue. They were first men- 
tioned in lAlcock et al.l (Il997al) and are commonly called "blue 
bumpers". 

We identify blue bumpers first from their position in the 
color-magnitude diagram and their low, chromatic amplification. 
For stars with (B - R)eros < 0.2 and i^eios < 18.6 (LMC) or 
T^eros < 18.9 (SMC), we require 



CIO: min(AR,AB) > 1-6 and 



A„-l 



1 



< 1.2 



where Ar and are the maximum amplifications observed in 
the two bands. These two amplifications are taken to be the fitted 
amplifications at the time of the maximum observed amplifica- 
tion. This value is used since the fitted amplification at f = fo 
may significantly overestimate a blue bumper amplitude if the 
peak region is not well sampled. 

There are 10 LMC light curves and 4 SMC light curves 
that pass cuts CI -CIO. Four of these light curves are most 
likely supernovae (SN) exploding in galaxies far behind the 
Magellanic Clouds. In the analysis of the F ull Sample of EROS- 
2 Magellanic light curves (Tisserand"2004), 31 such supernovae 
were found. All SN in our sample have fitted timescales in the 
range 25-50 days, asymmetric light curves with a faster rise time, 
and larger variations in the bluer passband. For a fraction of 



them, about 20%, the host galaxy is visible, which makes them 
indisputable SN. In order to identify the remaining SN, we have 
devised a fitting function with a time asymmetry parameter S ; 
the function reduces to simple microlensing for S - 0. The fit- 
ting function is obtained from simple microlensing ( fTOl l by re- 
placing the Einstein timescale Ie by the "supernova" varying 
timescale 



ts = tE 



1+5 arctan 



(13) 



For positive (negative) values of S , the rise time is faster (slower) 
than the decrease time. When fitting this function to the remain- 
ing light curves, we expect to find S > for the SN. Based on 
the study of the supernovae found in the Full Sample analysis 
and on simulated microlensing events, we choose to eliminate 
light curves with blue-band asymmetries by requiring 



Cll 



Sh < 0.30. 



This cut eliminates 3 light curves from the Bright-Star Sample, 
leaving 1 1 light curves, mostly with very small amplitudes, uq > 
1. Since the uq distribution for microlensing events is flat, the 
events with mo > 1 in Figure |6] must be mostly background. We 
therefore require 



C12 



Mo < 1.0 . 



This leaves two light curves. One of them is superimposed on a 
clear background galaxy. It is therefore most likely a supernova 
and we eliminate this event leaving one event to be discussed in 
the next section. 



5. Candidate sample and verifications 

After the cuts described in the previous two sections have been 
applied, only one candidate microlensing event remains in the 
Bright-Star Sample, EROS2-SMC-1. Its light curve is shown 
in fig. |7] and its characteristics are given on the first entry in 
Table IH It has been known sinc e 1997 (lAlcock et al.lll997bt 
iPalanque-Delabrouille et allll998i) and is identical to candidate 
MACHO-97-SMC-1. It has one of the two longest timescales 
(103 d) and the highest luminosity of all published microlensing 
candidates reported toward the Magellanic Clouds. The star is 
separated by 1.6arcs ec from another star that is 1.2 mag fainter 
(lUdalskiet al.lll997h causing ~ 30% blending in EROS-2 im- 
ages. Including the blend in the light curve fit increases fn to 
125d0The star displays a 4-5% variability with a period P - 
5.12 d jPal anque-Delabrouille et al. 19 98). For fur t her details on 
this event, see Palanque-Delabrouille et al.l (ll998l) . lAfonso et al.l 
(Il999ll2003al) and .Assef et aL(.2006) . 

As a check on the cuts CI - 12 leading to EROS2-SMC- 
1 and to search for nonstandard microlensing events (e.g. those 
due to binary lenses), a large number of light curves were visu- 
ally scanned. Among them were all events satisfying C1-C8 and 
Mo < 1, all events satisfying C1-C7 and mo < 0.5, all events sat- 
isfying CI -CIO, and all events satisfying CI -CI 2 but with R„ax 
increased by 0.2 mag. Only one interesting event was found, 
EROS2-LMC-15 shown in Figure [8] The beginning of the event 
is very similar to a standard microlensing event but subsequent 
points are too high[3 It was rejected by the X^,eak ^^^^ 



^ When blending is taken into account, this event has the longest 
of any published Magellanic event. 

OGLE-3 data confirm the high flux of the post-peak points (A. 
Udalski, private communication). 
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sm0054m 5761 EROS2-SMC-1 




Fig. 7. The light curves of EROS-2 microlensing candidate ER0S2-SMC-1 (star sm005-4m-5761). Also shown is the color- 
magnitude diagram of the star's CCD-quadrant and the excursion of the event. 




Fig. 8. The light curves of EROS-2 star kn057-0n-29305 (r.a.= 79.9488 deg, dec.= -70.7741 deg). Also shown is the color- 
magnitude diagram of the star's CCD-quadrant and the excursion of the event. 
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20.93 


0.23 


2204.1 


57.8 


0.36 


c Ar\ O /I £iC\ 

640.5/369 














1 Q 
ZZlj.J 


jl. 1 


U. 


ZyZ.L/ZZo 


EROS2-LMC-10 


lm045-5n-26323 


80.9925 


-65.8733 


20.55 


0.81 


1600.3 


37.1 


0.33 


580.0/505 














1602.2 


35.9 


0.51 


649.1/447 


ER0S2-LMC-11 


lm061-4m-15782 


88.4392 


-71.2536 


20.55 


0.55 


451.9 


51.3 


0.36 


456.5/444 














446.6 


46.1 


0.37 


458.7/394 


EROS2-LMC-12 


lm085-61- 14234 


89.8750 


-74.5676 


20.19 


0.34 


2664.0 


49.6 


0.27 


580.4/391 



Table 2. Events found in the analysis of the full sample of EROS-2 stars dTisserandlliOOj) . Only event EROS2-SMC-1 is in the 
Bright-Star Sample and passes the cuts C1-C12 and therefore is used for limits on the optical depth. All fits assume zero blending 
and no intrinsic variability. The two values of the fitted parameters to, t^ and mq are for Beros (first line) and /?eros (second line). (There 
are no /?ero.s points during the event for EROS2-LMC-12.) The time of maximum, to is given as JD-2450000. 



may, in fact, be a lensing event due to a binary lens. Note that its 
position in the color-magnitude diagram suggests that the star is 
either not in the LMC or behind a foreground Milky Way star. 

An important verification of the analysis of the Bright- 
Star Sample presented here comes from an analy sis of the 
Full S ample of Magellanic light curves described in iTisserandl 
(12004 . Because faint stars were analyzed, the cuts, given in 
AppendixlAl were generally slightly stricter than those described 
in SectionlD However, the analysis did not require observation of 
the event in both colors. It therefore provides an important check 
on our analysis. The candidates found in this analysis are listed 
in Table |2] and their light curves shown in Appendix [A] Two 
candidates EROS2-SMC-1 and 5 are in the Bright-Star Sample. 
Candidate 5 does not pass the cuts presented here because it fails 
C12. We also note that it shows non-microlensing-like variations 
in the light curve of the MACHO collaboratioiQ- 

Two other independent analyses were performed on the SMC 
data, restricted to the first five years of data. In both cases, the 
only event found in the Bright -Star Sample w as EROS2-SMC- 1 . 
The first analysis, reported in lAfonso et alj (l2003a ). was based 
on a larger set of stars (5.2 million); the number of analysis cuts 
common to it and the present analysis is small; the computer pro- 
grams were written independently. The second analysis followed 
a complete reprocessing using a new implementation of differen- 
tial photometry developed bv .Le Guillou (2Q03). The technique 
would allow us to find events not on cataloged stars. 

We attempted to check our efficiency for finding microlens- 
ing events by considering the eve nts published by the MACHO 
collaboration. The 13-17 events jAIcock et aljr2000bl) used by 
them to measure the optical depth toward the LMC are listed 
in Table [3] Only 2 of the 17 stars (MACHO-LMC- 1 8 and 25) 
are bright enough to be in our Bright-Star Sample. Three of 
the 17 events occurred after the beginning of EROS-2 opera- 
tions but none of these three events occurred on stars in our 
Bright-Star Sample. One of the three (MACHO-LMC- 15) was 
on a star too dim to be seen by EROS -20 The two other events 



" http://wwwmacho.mcmaster.ca/Data/MachoData.html. There is no 
information at this site on the EROS-2 LMC candidates of Table|2]since 
none are in MACHO fields. 

The 1.27 m telescope of the MACHO collaboration allowed them 
to use fainter stars than EROS-2 with its 1 m telescope. 



(MACHO-LMC- 14 and 20) are seen in the EROS-2 images and 
give microlensing parameters compatible with those measured 
by MACHO (see Table|3]l. 

MACHO-LMC- 14 was selected in the MACHO A analysis. 
In EROS-2, it was located in a defective zone of 200 x 2000 
pixels on CCD 7 of our blue camera, near the edge of the mo- 
saic. The corresponding star was not cataloged in this band, and 
the star failed the requirement to be observed in both passbands. 
Consequently, it does not appear in Table |2] 

MACHO-LMC-20 was selected in the MACHO B analysis. 
In EROS-2, it is located 35 arcsec from a very bright Galactic 
star (about V = 10 vs. 21 for the candidate). In this analysis, 
stars too near bright stars were eliminated so this candidate does 
not appear in Table |2] 



6. The status of published Magellanic microlensing 
candidates 

In this section, we review and update the status of the published 
Magellanic microlensing candidates of the EROS and MACHO 
collaborations. 

EROS has, in the past, used 1 1 candidates to place upper 
limits on the optical depth toward the LMC and to measure the 
depth toward the SMC. Six of the candidates are in the Bright- 
Star Sample considered here. The candidates and their present 
status are given in Table |4] 

Of the candidates in Table g] only EROS2-SMC-1 re- 
mains as a candidate microlensing event. The others have 
been eliminated either because continued observations of the 
same stars show further variability on the light curves, or be- 
cause improved photometry, in some cases complemented by 
spectroscopy, lead to re-interpreting the candidates as vari- 
able stars. Before the analysis presented here, the following 
EROS candidates were eliminated: EROSl-LMC-2, presented 
in A ubourg et al.l (Il993h and which displayed a new varia- 
tion 8 years later (iLasserre et al. 20()C ); and candidate EROS2- 
LMC-4, p resented in iLasserre et al.l ( 20001) and eliminated in 
iMilsztain et al. {200^. 

Before this analysis started, there were 5 surviving LMC 
microlensing candidates from EROS, one from EROS-1 (num- 
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MACHU candiclate 


bKUiS-i stai 










A 


comments 


MACHO-LMC- 1 


lm019-6k-14879 


18.97 


0.61 


-942.9 


17.20 


7.15 




MACHO-LMC-4 


lm057-01-27696 


19.85 


0.24 


-353.1 


22.70 


2.92 




MACHO-LMC-5 


lm057-0k-27469 


19.97 


1.52 


-976.1 


37.80 


47.28 




MACHO-LMC-6 


too faint 






-802.9 


45.80 


2.43 




MACHO-LMC-7 


lm010-2k-21383 


20.47 


0.52 


-536.5 


51.45 


5.91 




MACHO-LMC-8 


too faint 






-612.2 


33.20 


2.19 




ATA /^TT/~V T A T/^ O 1~» 1 ' 

MACHO-LMC-9 B-bmaty 


1 f\r\ 1 f\i /' a /" A 

ImOO 1-01-6864 


20.92 


0.83 


-400. 1 


89.60 


1.95 




ATA /^TT/~V T A T/^ t O 

MACHO-LMC- 1 3 


lm020- 6m- 19733 


21.03 


-0.07 


134.0 


50.05 


2.36 




ATA /~1TT/^ T A T/^ 1 A 

MACHO-LMC- 14 


1 r\r\'~\ '~i H /' 1 /' 

lm002-7n-24616 


18.84 




391.3 


50.05 (38.2) 


3.37 (2.77) 




ATA /~1TT/^ T A T/^ 1 C 

MACHO-LMC- 1 5 


too faint 






472.4 


18.40 


2.83 




ATA /^TT/~V T A T/^ t O 

MACHO-LMC- 1 8 


lm060-2n- 16057 


18.35 


0.47 


-217.2 


37.10 


1.54 


Bright Sample 


MACHO-LMC-20 B 


lm012-5n-23157 


20.09 


1.03 


llA.l 


36.35 (23.2) 


2.95 (3.1) 




MACHO-LMC-21 


dead zone 






-786.8 


46.60 


5.64 




MACHO-LMC-22 B 


dead zone 






-43.3 


114.65 


2.70 




MACHO-LMC-23 


lm055-3n-4994 


19.82 


0.37 


-237.7 


42.60 


2.41 


variable (Fig.llOl) 


MACHO-LMC-25 


lm017-2m-2056 


17.95 




-265.7 


42.60 


1.50 


Bright Sample 


MACHO-LMC-27 B 


lmOlO-On-23830 


19.00 


0.14 


-485.9 


25.25 


1.45 





Table 3. The 17 microlensing candidates of the MACHO collaboration jAlcock et al.ll2000bl) . Candidates 20, 22 and 27 (marked 
B) are low signal-to-noise candidates not satisfying the MACHO "A" requirements. Candidate 9 is due to a binary lens and does not 
satisfy their A requirements. The values of Ie and maximum amplification A,„fl^ are from the MACHO fit assuming no blending, and 
to is given as JD-2450000. Candidates 18 and 25 are each within 1 arcsec of an object of similar magnitude that is resolved by the 
MACHO analysis but blended in the EROS analysis. Candidates 14 and 20 occurred during EROS-2 observations, and the EROS-2 
values of ?£ and A^^^ are shown in parentheses (averages of /?eios and Beios measurements). 



Candidate 



EROS-2 star 



(B-R)„ 



original ref. 



status 



EROSl-LMC-1 lm058-2k-21915 18.75 0.34 

EROSl-LMC-2 lm043-6m- 9377 B-S 19.32 -0.04 

EROS2-LMC-3 lm034-61-20493 20.90 0.61 

EROS2-LMC-4 lm018-6n-23236 19.10 1.87 

EROS2-LMC-5 lm015-3n-22431 B-S 19.17 0.14 

EROS2-LMC-6 lm067-5m- 14700 21.01 0.63 

EROS2-LMC-7 lm070-3n-23389 21.00 0.76 

EROS2-SMC-1 sm005-4m-5761 B-S 18.13 -0.13 

EROS2-SMC-2 smOO 1-61- 13221 B-S 19.56 0.44 

EROS2-SMC-3 sm001-6n-16904 B-S 19.31 0.59 

EROS2-SMC-4 sm002-7m-21331 B-S 19.48 0.32 



Aubourg et al. (1993) 
Auboure et al. (1993) 
j.-asserre et al. (2000) 
Lasserre et al. (2000) 
Mlsztain et al. (2001) 
Milsztajn et al. (2001) 
Milsztain et al. (2001) 
Palanque-D. etal.(1998) 
Afonso et al. (2003a) 
Afonso et al. (2003a) 



2"'' variation (Tisserand 2004) (Figure |9]( 
2"'' variation (Lasserre et al. 2000) 
2"'' variation (Tisserand 2004) 
2'"' variation (Milsztain et a l. 2001) 
Supernova (Tisserand 2004) 
Supernova (Tisserand 2004) 
Supernova (Tisserand 2004) 
Fig.H] 

long period variable (Tisserand 2004) 
long period variable (Tisserand 2004) 
long period variable (Tisserand 2004) 



Afonso et al. (200 3a) 

Table 4. The 1 1 events of the EROS collaboration used in the past to set upper limits on the microlensing optical depth toward the 
LMC and to measure the depth toward the SMC. Candidates marked "B-S" in column 2 occur on stars in the EROS-2 Bright-Star 
Sample. All candidates except EROS2-SMC-1 have been eliminated as variable stars or as supernovae. 



ber EROS 1 -LMC- 1 ) and four from EROS-2 (numbered EROS2- 
LMC-3, 5, 6 and 7). 



The EROS-1 candidate EROS 1 -LMC- 1 displayed a new 
variation in the EROS-2 data (Figure|9]l in 1998, 6.3 years after 
the first one, of similar amplitude (a factor two) and timescale 
(27 days). This second variation is well fitted by a microlensing 
light curve. Because they are separated in time by more than 80 
Einstein timescales, the probability that these two bumps corre- 
spond to the microlensing of a double source star is lower than 
half a percent, even in the favorable case of the two stars being of 
equal luminosity. This candidate was thus rejected. Let us recall 
that it was already known to be a Be star (Beaulieu e t al. 199 5|) 
and w as thus already suspected of being variable dPaczviiskil 
I1996I) . Candidate EROS2-LMC-3 also displayed new variations 



between 1999 and 2002, of a more irregular nature, and was thus 
rejected. 

The hght curves of candidates EROS2-LMC-5, 6 and 7 
have been improved, due to the better template images in the 
pr esent analysis. The re duced photometric scatter, compared 
to iMilsztain et al.l (l200ll) . made apparent an asymmetry in rise 
and fall times; the asymmetry test (Cll) allowed us to iden- 
tify and reject them as supernovae. EROS2-LMC- 5 is identical 
to MAC HO-LMC-26, which had been rejected bv lAlcock et al.l 
(I2000bh for the same reason. 

The conclusion is that none of the former EROS LMC mi- 
crolensing candidates are still considered valid. Four displayed 
further variability, and three were identified as SN thanks to im- 
proved photometry. 
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There have been four EROS-2 SMC candidates, EROS2- 
SMC-1 discussed in Section |5l and EROS2-SMC-2, 3 and 4 
presented in Afonso et al. (2003a). Candidates 2, 3 and 4 were 
described as doubtful candidates, as all three display very long 
timescale variations, and look more like irregular variable stars. 
Another analysis of the sa me SMC data using differential pho- 
tometry (he Guillou"2003') allowed us to reduce the photomet- 
ric errors and subst antiate this in t erpreta tion. The additional two 
years of data since lAfonso et alj ( l2003ah have confirmed this, as 
further irregular variations have been observed. They are not se- 
lected in the present analysis, and are now considered as defini- 
tively rejected. 

We have also searched for further variations in the 17 
stars used by the MACHO collaboration to measure the opti- 
cal depth toward the LMC (Table O. One of the candidates, 
MAC HO-LMC-23, sh owed a further variation in the EROS-2 
data jTisserandl 12004'). 6.8 years after its first variation in the 
MACHO data. Its EROS-2 lightcurve is shown in Figure [TO] As 
such, we can eliminate it as a microlensing candidate. In spite 
of this, both the variation shown in Figure \W\ and the original 
variation in the MACHO data are quite achromatic, indicating 
that achromaticity is not a fool-proof crite rion for select i ng mi- 
crolensing events. We note however that iBennett et alj (l2005h 
argued that, even without considering its further variation seen 
in the EROS-2 data, the form of its light curve made MACHO- 
LMC-23 a weak microlensing candidate. 

The MACHO collaboration has also reported candidate 
events found by their alert system. The most notab le is 
MACH099-LMC-2 that was studied bv lBond et al] (|2002|) . This 
impressive microlensing candidate was on a star too faint to be 
in the EROS Full Sample. 

Besides the 17 LMC events of the MACHO collaboration, 
they have reported two candidates in the SMC, though they 
have reported no systemati c search for SMC events. The first 
MACHO SMC candidate ( Alcocket al.l ll997b) is identical to 
EROS2-SMC-1 (Section EJ. The second candidate, MACHO- 
98-SMC-l, concerned a star too faint to be included in our Full 
Sample of SMC stars. The event was detected by the alert sys- 
tem of the MACHO group in May 1998 and a probable caustic 
crossing due to a double lens was announced soon after. The 
source star was monitored by most microlensing groups, includ- 
ing EROS-2. The second caustic crossing was fully measured, 
which allowed the determination of the relative proper motion 
between the lens and source, allowi ng one to conclude t hat th e 
event was due to a lens in the SMC (lAfonso et alj|1998ll2000l) . 
H 

7. The detection efficiency 

To measure the optical depth from the detected events, or lim- 
its on this quantity, we first need to evaluate the detection effi- 
ciency as a function of the time scale f£. This was determined by 
using Monte Carlo simulated light curves : we superimpose ar- 
tificial microlensing events on a representative sample of light 
curves, coiTesponding to 2% of the Bright-Star Sample from 
each of the 98 monitored fields. The light curve for a simple mi- 
crolensing event (i.e. point-source point-lens zero-blending) is 

Note that there exists a third, unpublished, microlens- 
ing candidate toward the SMC, OGLE-2005-SMC-OOl 
(http://bulge.astro.princeton.edu/~ogle/). Its long timescale ( > 150 d), 
large amplification ( > 12) and bright source (/ = 18.2) offer the 
prospect of completely resolving the microlensing parameters degener- 
acy through measurements of lightcurve deformations, e.g. of that due 
to parallax. 



0.6 




0.0 1.0 2.0 3.0 



log(tE) 

Fig. 11. The detection efficiency for unblended (a = 1) mi- 
crolensing light curves ( fTOl i as a function of Ie for LMC (bold 
line) and SMC (light line) fields. The efficiency applies to a 
2500 d total observing period. The efficiency used for calcula- 
tion of the optical depth is the efficiency shown here, multiplied 
by 0.9 to take into account lensing by binary lenses. 
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Fig. 12. The ratio between the summed efficiency, e\(a\)+e2{a2) 
and the unblended efficiency e{a - 1) for pairs {ai,a2) taken 
from the artificial star analysis of Section [3] The ratio is a func- 
tion of fE. The four curves correspond to the four studied fields, 
lm009, 019, 034 and 048. 



described by three parameters of ( fTOl l: date of maximum amplifi- 
cation ?(), impact parameter mq and time scale fg. Blended events 
have the additional parameter a in (fTTT i. The microlensing pa- 
rameters are chosen at random: fo follows a flat distribution over 
our 2500 days observing period, JD 2,450,242 till 2,452,742; mq 
is picked randomly between and 1 .2; and tz is chosen at ran- 
dom from a distribution flat in ln(f£), between 1 and 1000 days. 
Each star in the 2% sample is actually used thrice in the simu- 
lation, once per decade in Ie- The simulation takes into account 
the relative variation of photometric errors. 

Simulated light curves were then fed into the analysis chain 
to find the fraction that are recovered by our detection algorithm. 
The detection efficiency in a given bin, eitz), is then given by 
the ratio between the number of events passing all selection cri- 
teria in this bin and the number of microlensing events generated 
in the same bin with mo < 1 • 

Figure [TT| shows the LMC and SMC efficiencies for un- 
blended events {a - 1) as a function of t^- For the range of 
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lm0582k 21915 EROSl-LMC-1 




1.0 



Fig. 9. The EROS-2 light curve of EROS- 1 microlensing candidate EROS 1 -LMC- 1 . The curve shows a second variation, 6.3 years 
after the variation observed in EROS-1. Also shown is the color-magnitude diagram of the star's CCD-quadrant and the excursion 
of the event. 



Im0553n 4994 MACHO-LMC-23 




Fig. 10. The EROS-2 light curve of MACHO microlensing candidate MACHO-LMC-23. The curve shows a second variation, 6.8 
years after the variation seen by MACHO. Also shown is the color-magnitude diagram of the star's CCD-quadrant and the excursion 
of the event. 
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Fig. 13. The detection efficiency for simple microlensing light 
curves (fTOl l as a function of mq (top) and of /?eros (bottom) for 
LMC and SMC events in the range 10 d < te < 200 d. 

Ie of interest for this work, the efficiency varies from ~ 0.25 at 
fE = 10 d to ~ 0.45 at Ie = 200 d. 

Blended events have an efficiency that is reduced because an 
impact parameter smaller than unity is necessary to produce an 
amplification of 1.34. The efficiency is also modified because 
the time scale is reduced. To sufficiently good approximation we 
find 



M,(2.18, a,) 
05 



e{a ^ 1, M,(l,a,)fE) 



(14) 



where Ui(A,„ax, ad is the impact parameter necessary to produce 
a maximum amplification A,„ax- 

To evaluate the summed detection efficiency for a realistic 
distribution of ai and 02 we used the pairs from the artificial 
images of Section[3] Figure|4] Figure[T2]shows the ratio between 
the calculated sum (for the four studied fields) and the unblended 
efficiency. For Ie = 40 d the efficiency is reduced by a factor 
ranging from 0.90 in the sparse fields lm048 and lm034 to 0.92 
and 0.97 in the denser fields lm019 and lm009. The brightest 
star, / = 1 accounts for 95% of the rate in the sparse fields and 
88% in the dense fields. 

Figure [T3] shows the efficiency for unblended events as a 
function of mq and of 7?eios for events with 10 d < tE < 200 d. 
The efficiency for the Bright-Star Sample has a much weaker 
dependence on uq and /?eros than that for the Full Sample. For 
the Full Sample the efficiency falls rapidly with increasing uq 
and /?eros (Tisserand 2004). 

The efficiencies in Figures [TT| and [13] are for the detection of 
microlensing events due to simple lenses. Events due to binary 




150 200 
te (days) 



250 



Fig. 14. The Ie distribution dT/dtE expected for IM© lenses 
in a spherically symmetric isoffiermal halo with core radius 
5kpc, i.e. ffie S m odel used by the MACHO collaboration 
(lAlcock etal.ll2000bl) . 



lenses with caustic crossings are discriminated against, mostly 
by C9. Of the 17 LMC events of the MACHO collaboration 
(lAlcock et al.ll2000 b) only 1 event, MACHO-LMC-9, is of this 
type and would not pass our selection criteria. This event cor- 
responds to ~ 10% of their optical depth. We note also that to- 
ward the Galactic Bulge ~ 10 % of ffie observed mi crolensing 
events are due to binary lenses (lUdalski et al.ll2000al) . To com- 
pensate for this loss of efficiency, we conservatively reduce the 
efficiency of Figure[TT]by a factor 0.9 when calculating limits on 
the optical depth. 



8. Limits on the abundance of machos 

The microlensing optical depth, r, is defined as the probability 
that any given star, at a given time, is amplified by at least 1.34, 
i.e. with an impact parameter u < I. From a set of A^ev events, r 
can be estimated from 



tE 



2A^starsrobs \ £ 



(15) 



where (fE/e) is the mean tE divided by efficiency for the observed 
events. 

In the LMC we have found no events so we can only give an 
upper limit on t by calculating the expected number of events 
as a function of t as given by For this an alysis we use the 
tE distribution of the S model ( Alcock et al.l HoOOb) shown in 
Figure [14] It relates {tE) to the macho mass (assumed unique): 
(tE) = lOdyfMjMQ. Limits using other halo models or macho 
mass distributions can be found to often good approximation by 
simply scaling (O with e({tE}) / {(e) ■ 

The expected number of LMC events for t/,„c = 4.7 x 10"^ 
as a function of lens mass, M, is shown in Figure [Tsh. For M = 
0.4Mo, we have (Je) = 44 d, (e) = 0.35, A^,,^^ = 5.5 x lO*" and 
Tohs = 2500 d, giving 32 LMC events for EROS-2. We add 7 
LMC events for EROS-1 to give a total of 39 expected events for 
n„, = 4.7 X 10-^. 

For no observed events (A^ < 3.0, 95% CL), the 95% CL 
upper limit on the optical depffi is 



(16) 



4.7 X 10-7 N,x(4.7 X 10-7) 
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For 39 expected events, The upper limit is then t;,,,^ < 0.36 x 
10"^. The limit on ti^c as a function of M is shown in Figure 
[TSb . In the /e range favored by the MACHO collaboration, we 
find 

Time < 0.36 X 10"^ X [1 +log(M/0.4Mo)] 95%CL , (17) 



150 



I.e. 



/ < 0.077 X [1 + log(M/0.4Mo)] 95%CL 



(18) 



where / = Tin]c/4.7 x 10"^ is the halo mass fraction within the 
framework of the S model. This limit on the optical depth is 
significantly below the value for the centr al region of the LMC 
measured by the MACHO collaboration (lAlcock et al.ll2000bl) . 
Timc/10~^ = 1.2;^Q3(sfaf.) ± 0.36(iyi.) and the revised value of 
lBennettl(l2005h .Ti^./10'^ = 1 .0 ± 0.3. The lAlcock et al.l (l2000bh 
optical depth used for the entire LMC predicts that EROS would 
see ~ 9 LMC events whereas none are seen. 

For the SMC, the one observed event corresponds to an opti- 
cal depth of 1.7 X 10"^ (A^star - 0.86 x 10*=). Taking into account 
only Poisson statistics on one event, 0.05 < A^obs < 4.74 (90% 
CL) this gives 



0.085 X 10" 



< 8.0 X 10" 



90%CL . 



(19) 



This is consistent with the expect ations of lensing by obj ects in 
the SMC itself, r,^^ ~ 0.4 x 10"^ dGraff & Gardinei]|l999,) . The 
value of fE - 125 d is also consistent with expectations for self- 
lensing {t^) ~ 100 d for a mean lens mass of O.35M0. 

We also note that the self-lensing interpretation is favored 
from the absence of an indication of parallax in the light curve 
dAssefet al.ll2006l) . 

We can combine the LMC data and the SMC data to give a 
limit on the halo contribution to the optical depth by supposing 
that the SMC optical depth is the sum of a halo contribution, 
Tsmc-haio = cTimc (tt ~ 1.4) and a self-lensing contribution Tsi. 
(We conservatively ignore contributions from LMC self-lensing 
and from lensing by stars in the disk of the Milky Way.) For one 
observed SMC event with Je = 125 d and zero observed LMC 
events, the likelihood function is 

^Tlnic,Tsl) oc [aTimcr[,(fE) + Tsir^l(fE)]eXp[-A^(Timc,Tsl)] 

where N{T\^i^,Ts\) is the total number of expected events (LMC 
and SMC) as a function of the two optical depths as calcu- 
lated with equation ([8]l. The function T'^(tE) is the distribu- 
tion (normalized to unit integral) expected for halo lenses of 
mass M (Figure [T4b and ^',.^{t ^ is the expected distribution for 
SMC self-lensing taken from iGraff & Gardineii (11999). We as- 
sume the SMC self-lens ing optical depth is that calculated by 
iGraff & Gardineii ( Il999h though the results are not sensitive to 
this assumption. For macho masses less than IM©, the likeli- 
hood function is maximized for ri^c = because there are 
no LMC events in spite of the greater number of LMC source 
stars. For M < O.IM© the limit on the halo contribution ap- 
proaches that one would calculate for no candidates in either 
the LMC or the SMC because the observed of 125 d is too 
long for a halo event. The calculated upper limit is shown as 
the dashed line in Figure [T5b. In the mass range favored by the 
MACHO collaboration, the limit is slightly lower than that us- 
ing only the LMC data. The combined limit would be somewhat 
stronger if we assumed an oblate halo (a < 1 .4) and somewhat 
weaker if we assumed a prolate halo (a > 1.4). Constraints on 
the shape of th e Milk y Way halo were recently summarized by 
iFellhauer et alj (l2006h who argued that the observed bifurcation 



.100 



II 



P 50 




2 2 

logM= 21og«tE>/70d) 



0.6[ 



0.4 



0.2 



0.0 





1 1 1 1 1 1 


1 




MACHO / 






95% cl / 






\ EROS-2 + EROS-1 / + 






\ upper limit (95% cl) 1 

1 1 1 1 i" 1 





-6 -4-2 2 

logM= 21og«tE>/70d) 

Fig. 15. The top panel shows the numbers of ex pected events 
as a fun ction of macho mass M for the S model of lAlcock et al.l 
(l2000bl) . The expectations for EROS-2 -LMC, SMC (this w ork) 
are shown along with those of EROS-1 (iRenault et alJl997h with 
contrib utions from the phot ographic plate progr am jAnsari et al.l 
Il996ah and CCD program (iRenault et al.lll998l) . The number of 
events for ER0S-2-SMC supposes Tsmc - 1.4T;„,f. In the lower 
panel the solid line shows the EROS 95% CL upper limit on 
/ = T"irnc/4.7 X 10"^ based on no observed events in the EROS- 
2 LMC data and the EROS-1 data. The dashed line shows the 
EROS upper limit on Time based on one observed SMC event in 
all EROS-2 and EROS-1 data assuming Tsme-haio - l-4Time. The 
MACH O 95% CL. curve is taken from Figure 12 (A, no Imc 
halo) of lAlcock et all (l2000bl) . 



of the Sagittarius Stream can be explained if the halo is close to 
spherical. 

A possible systematic error in our result could come from 
our assumption that the optical depth due to binary lenses is 
small, 10% of the total. An alternative strategy would have been 
to relax the cuts so as to include the event shown in Figure [8] 
We have chosen not to do this because the light curve itself is 
not sufficiently well sampled to establish the nature of the event 
(other than that it is not a simple microlensing event) and also 
because of its anomalous position in the color-magnitude dia- 
gram. We note also that the optical depth associated with the 
event, t = 0.7 x 10"**, is a factor ~ 4 below the upper limit (fTTl i. 

Another important question concerns the influence on our 
results of the Bright-Sample magnitude cut. Since the cut was 
not established before the event search, it is natural to ask if the 
position of the cut was chosen to give a strong limit. In fact, 
elimination of the cut would not change significantly the conclu- 
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sions of this paper. Four additional events (EROS2-LMC-8, 9, 10 
and 1 1 from Table |2]i were found with the analysis described in 
Appendix [a] The baselines of the source stars of these events 
are on average 0.9 magnitude below the Bright Sample cut. An 
additional event (EROS2-LMC-12) with no /?eros data during the 
variation was found. The values of fe for these events are in the 
range 10-60 days, similar to those in the MACHO sample. The 
optical depth associa ted wi th the four events seen in two colors 
is (1.7 ± 0.6) X 10"^ dTisser and 2004), not in contradiction with 
the limit from the Bright-Sample analysis. The 95% CL limit on 
the halo fraction is 12% at M = OAMq. We prefer to use the 
limit from the Bright Sample for the reasons already mentioned: 
better understanding of blending and a superior photometry al- 
lowing better rejection of variable stars. Indeed, inspection of the 
light curves of EROS2-LMC-8-12 in Figures lAJlfes] indicates 
that most are not especially convincing candidates. An excep- 
tion is EROS2-LMC-8 but its position in the color-magnitude 
diagram and strongly chromatic magnification indicate that the 
event is most likely due to a lens in the Milky Way disk. 

9. Discussion of tlie results 

The limits shown in Figure[T5]rule out machos as the majority of 
Galactic dark matter over the range 0.6 x IO^^Mq < M < 15Mq. 
The limits are / < 0.04 for IQ-^Mq <M < IO^'Mq and / < 0.1 
for 10"''Mo < M < IMq. We note that even stro nger limits in 
the range were found dAlcock et al.l 

11998 ) by combining the data from the EROS-1 CCD program 
with those of the MACHO program. These data gave / < 0.07 
for M ~ 3 X IO^^Mq. It should also be possible to improve limits 
in the range IOMq < M < IOOM0 by combining t he results pre- 
sented here with the MACHO high mass results (lAlcock et alJ 
I2001b^. This may narrow the small remaining allowed ma- 
cho mass range between the range excluded by microlensing 
and that excluded by the a bundance of halo wide binary stars 
jYoo. Chaname and Gouldll2004) . 

Initially, the EROS and MACHO programs were primar- 
ily motivated by the search for halo brown dwarfs of mass 
0.02 - O.O8M0. Such objects are clearly ruled out as primary 
components of the Milky Way halo by FigurefTSl (i mplying / < 
0.04) a nd the data of the MACHO collaboration (lAlcock et alJ 
I2000bl) . 

The observation of 17 events by MACHO with (fn) ~ 40d 
suggested the existence of machos of M ~ 0.4Mo. Such ob- 
jects could be invisible (e.g. primordial black holes) or faint (e.g. 
cool white dwarfs). However, the latter are not seen at the cor- 
responding level in multi-col or surveys fc ates et alll2004l) and 
high proper-motion surveys jGoldman etal.1 l2002l:lReidl I2005L 
and references therein). 

At any rate, the EROS limit ( fTTb is significantly less than 
the MACHO result ti^c = ^-^-ot ^ However, there are 
considerable differences between the EROS and MACHO data 
sets that may help to resolve the conflict. Generally speaking, 
MACHO uses faint stars in dense fields (1.1 x 10^ stars over 
13.4 deg^) while EROS-2 uses bright stars in sparse fields (0.7 x 
lO'' stars over Udeg^). Of these bright EROS-2 stai's, 0.2 x 10^ 
are in MACHO fields. 

The use of dense fields by the MACHO group suggests that 
the higher MACHO optical depth may be due, in part, to self- 
lensing in the inner parts of the L MC. This would contradict 
LMC models (M ancini et alJl200 4') which suggest that only 1- 
3 MACHO events should be expected to be due to self-lensing. 
In fact, MACHO-LM C-14 is known to be due to self-lensing 
(lAlcock etalJ2001cl) because it has a binary source and the form 



of the accompanying deformation of the lightcurve with respect 
to the simple microlensing lightcurve requires that the lens be in 
the LMC. A second event, MACHO-LMC-9, is due to a binary 
lens and the self-lensing interpretation can be avoided only by 
assuming that the source is also a binary system and that each of 
the two widely separated components happened to land on the 
caustic on the two successive observations made of the caustic 
entrance ( Alcock et al. 2000a). 

If it turns out that the self-lensing rate is higher than the 
model estimates in the MACHO fields but still negligible in the 
outer fields of the LMC, the disagreement between MACHO and 
this work is considerably reduced. Since only 1/3 of our Bright 
Sample stars are in MACHO fields, the EROS-2 95% CL up- 
per limit on t for the MACHO fi elds is the lim it (T% multiplied 
by a factor 3, consistent with the lAlcock et alJ fcOOOh) result as 
modified by Bennett (2001). 

A possible explanation for the discrepancy that is similar 
to self-lensing is the possibility that the events are due to halo 
lenses but the Halo is clumpy and that the MACHO-fields lie 
behind a clump of size less than that of the EROS-2 fields. The 
effect of a clumpy hal o on the optical depth was discussed by 
iHolopainen et al l (120061) though they did not discuss directly the 
possibility that it could resolve the EROS-MACHO controversy. 
At any rate, if this is the cause of the discrepancy, the EROS-2 
result gives the more representative optical depth because it is 
based on a larger solid angle. 

The use of faint stars by MACHO may also give an explana- 
tion of the disagreement. Only two of the 17 MACHO candidates 
(MACHO-LMC-18 and 25) are sufficiently bright to be in our 
Bright-Star Sample. Of these two, MACHO-LMC-18 is in the 
EROS Bright-Sample only because the EROS starfinder mixed 
two similar objects that are resolved by the MACHO starfinder. 

The use of faint stars by MACHO suggests two possible ex- 
planations for the disagreement. The first would be contamina- 
tion by variable stars that, in our Bright-Star Sample, are either 
not present or identified as such because of superior photomet- 
ric precision. The low photometric precision for the faint LMC 
stars makes most of the events less convincing than the events on 
bright stars in the Galactic Bulge or the one EROS event in the 
SMC. Indeed, one of the MACHO A events, MACHO-LMC- 
23, has been identified as a variable st ar by a second v ariation 
in the EROS- 2 data shown in Figure [TOl (lTisserandl2004l) . As al- 
ready noted. lBennett et all ^005) argued that MACHO-LMC-23 
was, in any case, a weak candidate and that its variability doesn't 
call into question the nature of the other MACHO candidates. 
Indeed, some of the MACHO candidates are very convincing 
microlensing candidates. In particular, MACHO-LMC 1, 5, 9, 14 
and 21 a re strong cand idates based only on MACHO photome- 
try while lBennettet al.l g.QQ5) argued that MACHO-LMC -4, 13 
and 15 are strong candidates because of high precision followup 
photometry. We might note, however, that of these candidates, 
MACHO-LMC-14 is most Ukely due to LMC self-lensing and 
MACHO-LMC-5 is due to l ensing by a normal red-d warf star in 
the disk of the Milky Way CAlcock et alJl2001dl:lGould 2004) . 

The second possible explanation related to the use of faint 
source stars is the large blending effects that must be understood. 
Alcock et al. (2000b, 2001a) suggest a 30% systematic error to 
reflect this uncertainty. The experience with the use of faint stars 
in the Galactic Bulge suggests that this uncertainty may be un- 
derestimated, though in principle the star distribution is better 
understood in the LMC than in the Bulge. At any rate, the exten- 
sion of the limits presented here, either by EROS or by OGLE-3 
or SuperMACHO, will require the use of faint stars and a good 
understanding of blending. 
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Appendix A: Candidate selection for the full sample 

In this paper we have concentrated on the analysis of the Bright- 
Star Sample of stars because this leads to the most reliable limits 
of the optical depth. A search for micro lensing events o n the 
full sample of stars was also performed dTisserandl |2004|) . The 
candidates found in this search are listed in Table|2]and displayed 
in Figures lA.lllA.SI In this section we list the selection criteria 
leading to this set of candidates. 

The criteria are very similar to those described in Section |4] 
for the Bright-Star Sample. However, the problems encountered 
with analyzing low precision light curves necessarily led to more 
complicated criteria to avoid the many spurious events caused by 
photometric problems. In all, 17 criteria, cl - cl7, were applied 
compared to the 12 criteria, Cl - C12, applied to the Bright-Star 
Sample. 

The first four criteria are identical to those applied to the 
Bright-Star Sample: 

cl = Cl c2 = C2 c3 = C3 c4 = C4 . 

We then required that the total amplitude of luminosity vari- 
ation along the light curve be greater than 5 times the point-to- 
point dispersion, crease in the light curve, recomputed after ex- 
cluding the most significant excursions: 

Here Fsmax and F5„„„ are the maximum and minimum fluxes av- 
eraged over any 5 neighboring measurements. 

For excursions with a regular variation, the point-to-point 
dispersion of the measured fluxes within the excursion, cr^„p,e.ic, 
obtained from the comparison of each measured flux with the 
linear interpolation of its two neighbors (in time), is normally 
smaller than the global dispersion, crj. j(,, i.e. the width of the dis- 
tribution of all fluxes within the excursion. To exclude irregular 
variations, we require that their ratio be sufficiently small; 

c6 : < 0.90 . 

O-exc 

Cuts c4 - c6 are applied independently in the two passbands. 
Light curves are retained if they are selected in one passband 
at least. This reduces the star sample to slightly less than 0.1% 
(28,500 objects) of the full sample. 

The next two criteria are similar to C5 and C9 applied to the 
Bright-Star Sample 



cl 



c8 : 




< 15 



< 10, 



peak 



changed the PSF. Faint stars near a diffraction feature in the PSF 
of bright stars were affected; this was seen as light curves with 
two plateaus, one (higher) before and one (lower) after May 
1998. A rather complicated algorithm was developed to iden- 
tify such light curves. It relies mostly on the relative flux val- 
ues of the two stars, their distance and the relative height of 
the two plateaus in the light curve. The details can be found in 
(Tisserand 2004.) but we note that it affects our efficiency only 
for/E > 150 d. 



c9 



no influence of telescope realignment 



At this point we accepted only light curves that passed crite- 
ria cl - c9 in both colors: 

clO : cl - c9 red and blue 

Criterion cl 1 is similar to C7: 

cl 1 : Npeak > 4 ; > 2 points in rise ; > 2 points in fall . 

The remaining cuts follow closely the cuts for the Bright- 
Star Sample 



cl2 = C6 



cl3 



cl4 



UQr < 1 or Uqi, < 1 

1 



2 2 



> 50 inRoi B 



The next cut, c9, deals with a background connected to the 
realignment of the telescope optics in May 1998, which slightly 



Criterion cl5 eliminates Ught echos from SN1987a by requiring 

cl5 : Distance event - SN1987a > 30 arcmin . 

Finally, cl6 and cl7 are the blue bumper and supernovae crite- 
ria CIO and Cll supplemented by a visual rejection of events 
superimposed on background galaxies. 

The first six events in Table |2] passed the cuts cl-cl7. A 
monochromatic analysis was performed where c 10 was removed 
and candidates were only required to pass the other criteria in ei- 
ther of the two colors. About 500 events were scanned by eye. 
Only one event, EROS2-LMC-12, was found, the others being 
due to photometric problems or long period variable stars. 
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Fig. A.l. The light curves of EROS-2 star lm055-7m-23303. Also shown is the color-magnitude diagram of the star's CCD-quadrant 

and the excursion of the event. 
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Fig. A.4. The light curves of EROS-2 star lm061-4m-15782. Also shown is the color-magnitude diagram of the star's CCD-quadrant 
and the excursion of the event. 
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Fig. A.5. The light curves of EROS-2 star lm085-61- 14234. Also shown is the color-magnitude diagram of the star's CCD-quadrant 
and the position of the source's baseline. 



Ansari, R.. Cavalier, R, Moniez, M., et al. (EROS-1 coll.) 1996a, A&A, 314, 94 

Ansari, R. et al. (EROS-2 coll.) 1996b, Vistas in Astronomy, 40, 519 

Assef, R.J., Gould, A., Afonso, C, et al 2006, ApJ, 649, 954 

Aubourg, E. et al., (EROS-1 coll.), 1993, Nature, 365, 623 

Bauer, F. et al. (EROS-2 coll.) 1997, Proceeding of the "Optical Detectors for 

Astronomy" workshop, ESO 
Beaulieu, J.-Ph., Ferlet, R., Orison, Ph. et al., (EROS-1 coU.), 1995, A&A, 299, 

168 

Bennett, D. P, Becker, A. C, & Tomaney, A. 2005, ApJ, 631, 301 
Bennett, D. P 2005, ApJ, 633, 906 
Bissantz, N. and Gerhard, O., 2002, MNRAS, 330, 591 
Bond, I.A., Rattenbury, N.J., Skuljan, J., et al. 2002, MNRAS, 333, 71-83. 
Calchi Novati, S., Paulin-Hemiksson, S., An, J., et al. (AGAPE coll.) 2005, 
A&A, 443, 911 

de Jong, J. T. A., Widrow, L.M., Cseresnjes, P, et al. (MEGA coU.) 2006, A&A, 
446, 855 

Evans, N.W. & Belokurov, V. 2002, ApJ, 567, LI 19 

Fellhauer, M., Belokurov, V., Evans, N.W., et al. 2006, ApJ, 651, 167 

Gates, E., Gyuk, G., Harris, H. C. et al. (SDSS coU.) 2004, ApJ, 612, 132 

Goldman, B., Afonso, C, Alard, C, et al. (EROS-2 coll.) 2002, A&A, 389, 69 

Gould, A. 2004, ApJ, 606, 319 

Graff, D. & Gardiner, L.T. 1999, MNRAS, 307, 577 

Griest, K. 1991, ApJ, 366, 412 

Hamadache, C, Le Guillou, L., Tisserand, R, et al (EROS-2 coll.) 2006, A&A, 
454, 185 

Han, C. and Gould, A. 2003, ApJ. 592, 172 

Holopainen, J., Flynn, C, Knebe, A., et al., 2006, MNRAS, 368, 1209 

Joshi, Y.C., Pandey, A.K., Narasimha, D. & Sagar, R. (Namital coll.) 2005, 

A&A, 433, 787 

Lasserre, T. 2000, Thesis Universite de Paris VI. http://tel.ccsd.cnrs.fr/ 
Lasserre, T, Afonso, C, Albert, J.N.et al., (EROS-2 coll.), 2000, A&A, 355, L39 
Le Guillou, L. 2003, Thesis Universite de Paris VI, http://tel.ccsd.cnrs.fr/ 
Mancini, L., Calchi Novati, S., Jetzer. Ph., & Scarpetta, G. 2004, A&A, 427, 61 
Milsztajn, A. and Lasserre, T. (for the EROS-2 coll.), 2001, Nucl. Phys. B (Proc. 

Suppl)91,413 
Paczyhski, B. 1986, ApJ, 304, 1 
Paczyiiski, B., 1996, ARA&A, 34, 419 



Palanque-Delabrouille, N., Afonso, C, Albert, J. N., et al. (EROS-2 coU.), 1998, 
A&A, 332, 1 

Petrou, M., 1981. Ph.D. thesis. University of Cambridge 

Popowski, P, Griest, K., Thomas, C. L., et al. (MACHO coll.) 2005, ApJ, 631, 
879 

Reid, I.N. 2005, ARA&A, 43, 247 

Renault, C, Afonso, C. Aubourg, E., et al. (EROS-1 coh.) 1997, A&A, 324, 
L69 

Renauh, C, Aubourg, E., Bareyre, R, et al. (EROS-1 coll.) 1998, A&A, 329, 522 
Riffeser, A., Fliri, J.. Bender, R., et al. (WeCapp coll.) 2003, ApJ, 599, L17 
Robin, A., Reyle, C, Demere, S., & Picaud, S. 2003, A&A, 409, 523; and 2004, 

A&A, 416, 157 
Sackett, PD. and Gould, A., 1993, ApJ, 419, 648 

Sumi, T., Wozniak, R R., Udalski, A., et al. (OGLE coll.) 2006, ApJ, 636, 240. 
Tisserand, P. 2004, Thesis Universite de Nice, http://tel.ccsd.cnrs.fr/ 
Udalski, A., Szymahsi, M., Kumiak, M., et al. (OGLE coll.), 1997, Acta Astron., 
47, 431 

Udalski A., Zebrun K., Szymanski M., et al. (OGLE coll.) 2000a, Acta Astron., 
50, 1 

Udalski, A., Szymanski, M., Kubiak, M., et al. (OGLE coll.) 2000b. Acta 
Astron., 50, 307 

Uglesich, R.R., Crotts, A. P S., Baltz, E.A., et al. (VATT coU.) 2004, ApJ, 612, 
877 

Wood, A. and Mao, S. 2005, MNRAS, 362, 945 

Yoo, J., Chaname J., and Gould, A. 2004, ApJ, 601, 311 

Zaritsky Z., Harris J., Thompson I., and Grebel, E. 2004, AJ, 128, 1606 



